Shot Noise in Graphene 



L. DiCarlo, 1 J. R. Williams, 2 Yiming Zhang, 1 D. T. McClure, 1 and C. M. Marcus 1 

1 Department of Physics, Harvard University, Cambridge, MA 02138, USA 
2 School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA 

(Dated: February 21, 2008) 

We report measurements of current noise in single- and multi-layer graphene devices. In four 
single-layer devices, including a p-n junction, the Fano factor remains constant to within ±10% 
upon varying carrier type and density, and averages between 0.35 and 0.38. The Fano factor in a 
multi-layer device is found to decrease from a maximal value of 0.33 at the charge-neutrality point 
to 0.25 at high carrier density. These results are compared to theories for shot noise in ballistic and 
disordered graphene. 



Shot noise, the temporal fluctuation of electric current 
out of equilibrium, originates from the partial transmis- 
sion of quantized charge [I]. Mechanisms that can lead 
to shot noise in mcsoscopic conductors include tunneling, 
quantum interference, and scattering from impurities and 
lattice defects. Shot noise yields information about trans- 
mission that is not available from the dc current alone. 

In graphene 13], a zero-gap two-dimensional semi- 
metal in which carrier type and density can be con- 
trolled by gate voltages |4j , density-dependent shot-noise 
signatures under various conditions have been investi- 
gated theoretically (5) [5]. For wide samples of bal- 
listic graphene (width-to-length ratio W/L > 4) the 
Fano factor, T ', i. e., the current noise normalized to the 
noise of Poissonian transmission statistics, is predicted 
to be 1/3 at the charge-neutrality point and ~ 0.12 
in both electron (n) and hole (p) regimes [5]. The 
value T = 1 — l/\/2 ss 0.29 is predicted for shot noise 
across a ballistic p-n junction [3]. For strong, smooth 
"charge-puddle" disorder, theory predicts T rs 0.30 both 
at and away from the charge-neutrality point, for all 
W/L > 1 7 . Disorder may thus have a similar ef- 
fect on noise in graphene as in diffusive metals, where 
T is universally 1/3 M M EEU III HI regardless of 
shape and carrier density. Recent theory investigates 
numerically the evolution from a density-dependent to 
a density-independent T with increasing disorder [14] . 
To our knowledge, experimental data for shot noise in 
graphene has not yet been reported. 

This Letter presents an experimental study of shot 
noise in graphene at low temperatures and zero mag- 
netic field. Data for five devices, including a locally- 
gated p-n junction, are presented. For three globally- 
gated, single-layer samples, we find T ~ 0.35 — 0.37 in 
both electron and hole doping regions, with essentially 
no dependence on electronic sheet density, n s , in the 
range \n 8 \ < 10 12 cm -2 . Similar values are obtained for 
a locally-gated single-layer p-n junction in both unipo- 
lar (n-n or p-p) and bipolar {p-n or n-p) regimes. In a 
multi-layer sample, the observed T evolves from 0.33 at 
the charge-neutrality point to 0.25 at n s ~ 6 x 10 12 cm~ 2 . 

Devices were fabricated by mechanical exfoliation of 
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FIG. 1: (a) Differential resistance R of sample Al as a function of 
back-gate voltage Vb g at electron temperature T e = 0.3 K, perpen- 
dicular field B± = 0, and source-drain voltage V s d = 0. (b) Dif- 
ferential two-terminal conductance g(V s( j = 0) as a function of B± 
and Vb g in the quantum Hall regime, after subtracting a quadratic 
fit at each B±. Lines of constant filling factors 6, 10, 14, and 18 
(dashed lines) indicate a single-layer sample, (c) Equivalent cir- 
cuit near 1.5 MHz of the system measuring current noise using 
cross correlation of two channels |16| . Current bias I contains a 
7.5 nA rma , 20 Hz part for lock-in measurements and a controllable 
dc part generating the dc component of V a ^ via the shunt resistance 
r = 5 kO. False-color scanning electron micrograph of a three-lead 
pattern defining two devices similar to Al and A2. Purple indicates 
single-layer graphene and gold indicates metallic contacts. 

highly-oriented pyrolytic graphite [4 . Exfoliated sheets 
were deposited on a degenerately-doped Si substrate 
capped with 300 nm of thermally grown SiC>2. Re- 
gions identified by optical microscopy as potential single- 
layer graphene were contacted with thermally evapo- 
rated Ti/Au leads (5/40 nm) patterned by electron-beam 
lithography. Additional steps in the fabrication of the p-n 
junction device are detailed in Ref. [15]. Devices were 
measured in two 3 He cryostats, one allowing dc (lock-in) 
transport measurements in fields |-B_i_| < 8 T perpendicu- 
lar to the graphene plane, and another allowing simulta- 
neous measurements of dc transport and noise |16j near 
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1.5 MHz, but limited to B± ~ 0. 

Differential resistance R = dV s d/dI (I is the current, 
and V s d is the source-drain voltage) of a wide, short sam- 
ple [Al, (W, L) — (2.0, 0.35) /mi] is shown as a function of 
back-gate voltage Vt, g at V^d = and B± = in Fig. 1(a). 
While the width of the peak is consistent with Al being 
single-layer graphene [TTJ, [TH] , more direct evidence is ob- 
tained from the QH signature shown in Fig. 1(b). The 
grayscale image shows differential conductance g = l/R 
as a function of Vb g and B±, following subtraction of 
the best-fit quadratic polynomial to <?(Vt, g ) at each B± 
setting to maximize contrast. Dashed lines correspond 
to filling factors n s h/eB± = 6, 10, 14, and 18, with 
n s = a(Vb g +l.l V) and lever arm a = 6.7xl0 10 cm~ 2 /V. 
Their alignment with local minima in 5g{V\, g ) identifies 
Al as single-layer graphene [HO |2D]. The Drude mean 
free path I = h/2e 2 ■ a /hp [2"T] . where fcp = -y/7r|n s |, is 
found to be ~ 40 nm away from the charge-neutrality 
point using the B± = conductivity a — (RW/L)^ 1 
[Fig. 2(a) inset]. 

Current noise spectral density Si is measured using 
a cross-correlation technique described in Ref. [TB] [see 
Fig. 1(c)]. Following calibration of amplifier gains and 
electron temperature T e using Johnson noise thermom- 
etry (JNT) for each cooldown, the excess noise Sf = 
Si — ik B T e g(V s< i) is extracted. (T4d) f° r sample Al is 
shown in Fig. 2(a). Linearity of Sj at high bias indicates 
negligible extrinsic (1// or telegraph) resistance fluctu- 
ations within the measurement bandwidth. For these 
data, a single-parameter fit to the scattering-theory form 
(for energy-independent transmission) |221 123) . 



SJ = 2eIT 



coth 



eV sd 
2k B T e 



2k B T e 
eV^d 



(1) 



gives a best-fit Fano factor T — 0.349. Simultaneously 
measured conductance g ~ 22.2 e 2 /h was independent 
of bias within ±0.5% (not shown) in the \V s d\ < 350 /iV 
range used for the fit. Note that the observed quadratic- 
to-linear crossover agrees well with that in the curve fit, 
indicating weak inelastic scattering in Al jTTJ [IS], and 
negligible series resistance (e.g., from contacts), which 
would broaden the crossover by reducing the effective V'sd 
across the sample. 

Figure 2(b) shows similarly measured values for T as 
a function of Vb g . T is observed to remain nearly con- 
stant for \n s \ < 10 12 cm~ 2 . Over this density range, 
the average T is 0.35 with standard deviation 0.01. The 
estimated error in the best-fit T at each Vbg setting is 
±0.002, comparable to the marker size and smaller than 
the variation in T near Vb s — 0, which we believe re- 
sults from mesoscopic fluctuations of T . Nearly identical 
noise results (not shown) were found for a similar sample 
(B), with dimensions (2.0,0.3) /im and a QH signature 
consistent with a single layer. 

Transport and noise data for a more square single- 
layer sample [A2, patterned on the same graphene sheet 
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FIG. 2: (a) Inset: Conductivity a = (RW/ L) — 1 calculated using 
i?(Vb g ) data in Fig. 1(a) and W/L = 5.7. Solid black circles cor- 
respond to o"(V s( j = 0) at the Vb g settings ol noise measurements 
shown in (b). Main: Excess noise 5| as function of V a( j near the 
charge-neutrality point, Vb g = —0.75 V. The solid red curve is the 
single-parameter best fit to Eq. (1), giving Fano factor T = 0.349 
(using T e = 303 mK as calibrated by JNT). (b) Best-fit T at 25 
Vbg settings across the charge-neutrality point for electron and hole 
densities reaching |n s | ~ 1.4 X 10 12 cm -2 , (c) R (left axis) and 
<t (right axis) of sample A2 as a function of Vbg (W/L = 1.4), 
with V s d = 0, at 0.3 K (solid markers) and at 1.1 K (open mark- 
ers), (d), (e) Crossover width T w (normalized to JNT-calibrated 
T e ) and T. obtained from best-fits using Eq. (1) to S'KVgd) data 
over |V sd | < 350(650) for T e = 0.3(1.1) K. 



as Al, with dimensions (1.8, 1.3) /jm] at T e = 0.3 K 
(solid circles) and T e — 1.1 K (open circles) are shown 
in Figs. 2(c-e). At both temperatures, the conductivity 
shows cr m i n a 1.5 e 2 /h and gives £ ~ 25 nm away from 
the charge-neutrality point. That these two values dif- 
fer from those in sample Al is particularly notable as 
samples Al and A2 were patterned on the same piece of 
graphene. Results of fitting Eq. (1) to Sj(V s d) for sample 
A2 are shown in Figs. 2(d) and 2(e). To allow for pos- 
sible broadening of the quadratic-to-linear crossover by 



3 



series resistance and/or inelastic scattering, we treat elec- 
tron temperature as a second fit parameter (along with 
!F) and compare the best-fit value, T w , with the T e ob- 
tained from Johnson noise. Figure 2(d) shows T w track- 
ing the calibrated T e at both temperatures. Small devia- 
tion of T w /T e from unity near the charge-neutrality point 
at T e — 0.3 K can be attributed to conductance varia- 
tions up to ±20% in the fit range |V s d| < 350 fiV at these 
values of Vbg- As in sample Al, T is found to be indepen- 
dent of carrier type and density over \n s \ < 10 12 cm , 
averaging 0.37(0.36) with standard deviation 0.02(0.02) 
at T e — 0.3(1.1) K. Evidently, despite its different aspect 
ratio, A2 exhibits a noise signature similar to that of Al. 

The lack of independence in T suggests that bias- 
dependent electron heating in the metallic reservoirs [T5] 
is negligible for our samples. This heating, originat- 
ing from imperfect dissipation of the generated power 
V 2 d /R, can affect shot noise measurements since these 
require |T4d| several times the thermal voltage (here, 
e\V s d\/kBT e < 10). In the presence of heating, fitting the 
excess noise Si(V sd ,T e + AT e (V sd )) ~ 4k B T e g to Eq. (1) 
overestimates T . The nearly equal values of T observed 
in Al and A2 despite the factor ^10 difference in R at 
comparable n s suggest that heating in the reservoirs is 
negligible [23] . 

Transport and noise measurements for a single-layer 
graphene p-n junction [15], sample C, are shown in Fig. 3. 
The color image in Fig. 3(a) shows differential resistance 
R as a function of Vbg and local top-gate voltage V tg . 
The two gates allow independent control of charge densi- 
ties in adjacent regions of the device [see Fig. 3(c) inset]. 
In the bipolar regime, the best-fit T shows little den- 
sity dependence and averages 0.38, equal to the average 
value deep in the unipolar regime, and similar to results 
for the back-gate-only single-layer samples (Al, A2 and 
B). Close to charge neutrality in either region (though 
particularly in the region under the top gate), Sj(V s d) 
deviates from the form of Eq. (1) (data not shown). This 
is presumably due to resistance fluctuation near charge 
neutrality, probably due mostly to mobile traps in the 
AI2O3 insulator beneath the top gate. 

Measurements at 0.3 K and at 1.1 K for sample D, of 
dimensions (1.8, 1.0) /im, are shown in Fig. 4. A ~ 3 nm 
step height between Si02 and carbon surfaces measured 
by atomic force microscopy prior to electron-beam lithog- 
raphy [25 suggests this device is likely multi-layer. Fur- 
ther indications include the broad R(Vb s ) peak [26] and 
the large minimum conductivity, er m i n ~ 8 e 2 /h at 
B± = [Fig. 4(a)], as well as the absence of QH sig- 
nature for \Bj_\ < 8 T at 250 mK (not shown). Two- 
parameter fits of Sj(14d) data to Eq. (1) show three no- 
table differences from results in the single-layer samples 
[Figs. 4(b) and 4(c)]: First, T shows a measurable de- 
pendence on back-gate voltage, decreasing from 0.33 at 
the charge-neutrality point to 0.25 at n s ~ 6 x 10 12 cm~ 2 
for T e = 0.3 K; Second, T decreases with increasing tem- 
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FIG. 3: (a) Differential resistance R of sample C, a single-layer 
p-n junction, as a function of back-gate voltage Vbg and top-gate 
voltage Vtg ■ The skewed-cross pattern defines quadrants of n and p 
carriers in regions 1 and 2. Red lines indicate charge-neutrality lines 
in region 1 (dotted) and region 2 (dashed), (b) 5J(V s( j) measured in 
n-p regime with (Vbgi Vtg) = (5, —4) V (solid dots) and best fit to 
Eq. (1) (red curve), with T = 0.36. (c) Main: Best-fit T along the 
cuts shown in (a), at which n s \ ~ n s 2 (purple) and n s i ~ —4 n s 2 
(black). Inset: Schematic of the device. The top gate covers region 
2 and one of the contacts. 



perature; Finally, T w /T e is 1.3-1.6 instead of very close 
to 1. We interpret the last two differences, as well as the 
sublinear dependence of S] on V s d (see Fig. 4 inset) as 
indicating sizable inelastic scattering [HH] in sample D. 
(An alternative explanation in terms of series resistance 
would require it to be density, bias, and temperature de- 
pendent, which is inconsistent with the independence of 
g on V sd and T e ). 

Summarizing the experimental results, we find that in 
four single-layer samples, T is insensitive to carrier type 
and density, temperature, aspect ratio, and the presence 
of a p-n junction. In one multi-layer sample, T does 
depend on density and temperature, and S]{V a d) shows a 
broadened quadratic-to-linear crossover and is sublinear 
in Vsd at high bias. We may now compare these results to 
expectations based on theoretical and numerical results 
for ballistic and disordered graphene. 

Theory for ballistic single-layer graphene with 
W/L > 4 gives a universal T — 1/3 at the charge- 
neutrality point, where transmission is evanescent, and 
T ~ 0.12 for \n s \ > 7r/L 2 , where propagating modes 
dominate transmission [5] . While the measured T at the 
charge-neutrality point in samples Al and B (W/L = 5.7 
and 6.7, respectively) is consistent with this prediction, 
the absence of density dependence is not: ir/L 2 ~ 3 x 
10 9 cm~ 2 is well within the range of carrier densities 
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FIG. 4: (color) (a) Differential resistance R (left axis) and conduc- 
tivity a (right axis) of sample D as a function of Vbg, with V* s( j = 0, 
at 0.3 K (solid markers) and at 1.1 K (open markers). (b),(c) Best- 
fit T w (normalized to JNT-calibrated T e ) and T to S^Vad) data 
over jy sd | < 0.5(1) mV for T e = 0.3(1.1) K. Inset: Sublinear de- 
pendence of Sj on V s d is evident in data taken over a larger bias 
range. Solid red curve is the two-parameter best fit of Eq. (1) over 
|V sd | < 0.5 mV. 



covered in the measurements. Theory |27] for ballistic 
graphene contacted with finite-density leads finds slight 
increments of T from 1/3 at the charge-neutrality point, 
in agreement with this experiment. However, T in this 
contact model remains density dependent. Theory for 
ballistic graphene p-n junctions [5] predicts T fts 0.29, 
lower than the value ~ 0.38 observed in sample C in 
both p-n and n-p regimes. We speculate that these dis- 
crepancies likely arise from the presence of disorder. Nu- 
merical results for strong, smooth disorder [7] predict a 
constant T at and away from the charge-neutrality point 
for W/L > 1, consistent with experiment. However, the 
predicted value T « 0.30 is ~ 20% lower than observed 
in all single-layer devices. Recent numerical simulations 
P3] of small samples (L = W ~ 10 nm) investigate the 
vanishing of carrier dependence in T with increasing dis- 
order strength. In the regime where disorder makes T 
density-independent, the value T ~ 0.35 — 0.40 is found 
to depend weakly on disorder strength and sample size. 

Since theory for an arbitrary number of layers is not 
available for comparison to noise results in the multi- 



layer sample D, we compare only to existing theory for 
ballistic bi-layer graphene [28] . It predicts T = 1/3 over 
a much narrower density range than for the single layer, 
and abrupt features in T at finite density due to trans- 
mission resonances. A noise theory beyond the bi-layer 
ballistic regime may thus be necessary to explain the ob- 
served smooth decrease of T with increasing density in 
sample D. 
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